Hydrodynamic Turbulence Has Infinitely Many Anomalous Dynamical Exponents 
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On the basis of the Navier-Stokes equations we develop the 
statistical theory of many space-time correlation functions of 
velocity differences. Their time dependence is not scale in- 
variant: n-order correlations functions exhibit n — 1 distinct 
decorrelation times that are characterized by n — 1 anomalous 
dynamical scaling exponents. We derive exact scaling rela- 
tions that bridge all these dynamical exponents to the static 
anomalous exponents £ n of the standard structure functions. 
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Experimental investigations of the statistical objects 
that characterize the small scale structure of turbulent 
flows are almost invariably based on a single point mea- 
surement of the velocity field as a function of time [Q-Q ■ 
The Taylor "frozen turbulence" hypothesis is then used 
to surrogate time for space. The results of this type of 
measurements are "simultaneous" correlation of the ve- 
locity field itself or of velocity differences across a scale R 
(structure functions) , or of velocity gradient fields like the 
dissipation field. Theoretical analysis which starts with 
the Navier-Stokes equations, on the other hand, states 
unequivocally that a closed form theory for the simulta- 
neous many space-point correlation functions of velocity 
differences is not available. The perturbative theory of 
turbulence is very clear about this: attempting to de- 
rive equations for simultaneous correlation functions one 
finds integrals over time differences of space-time corre- 
lation and response functions j^]. There are no small 
parameters like a ratio of time scales (as in turbulent ad- 
vection || ) or a small interaction (like in weak turbulence 
0) that allow a reduction of such a theory to a closed 
scheme in terms of simultaneous objects only. The aim of 
this Letter is to initiate a nonperturbative analytic theory 
of n-order space-time correlation functions, and to find 
the characteristics of such a theory that can be related 
to objects that are known from standard experiments. It 
should be pointed out that in addition to the fundamen- 
tal interest of such a theory it has important applications 
for the theory of scalar advection ; such a theory is 
beyond the scope of this Letter and is only mentioned as 
an additional motivation for the present analysis. 

In recent theoretical work 10 about the same time 
statistics of turbulence it was made clear why it is ex- 
tremely advantageous to consider many point correla- 
tion functions. The basic field under study is the differ- 
ence of the Eulerian velocity field u(r, t) across a scale 
R = r' r: w(r, r', t) = u(r', t) — u(r, t). The funda- 
mental statistical quantities are the simultaneous "fully 



unfused" n-rank tensor correlation function of velocity 
differences: 



F n (n,r' 1 ;r 2 ,r' 2 ; . . -;r n ,r' n ) 
= (w{r ll r' ll t)w{r 2l r 2l t) . . .w(r n ,r' n ,t)) 



(1) 



where pointed brackets denote the cumulant part of the 
average over a (time-stationary) ensemble. In this quan- 
tity all the coordinates are distinct. The more usual 
structure function S n (R) 



S n (R) = (\ w(r,r',t)\ n ) 



R = r' r 



(2) 



is obtained by fusing all the coordinates 7% into one point 
r, and all the coordinates r[ into another point r + R. 
Obviously in using the functions of many variables F n 
instead of the one variable function S n (R) one is paying 
a heavy price. On the other hand this has an enormous 
advantage: when one develops the theory for S n {R) on 
the basis of the Navier-Stokes equations one encounters 
the notorious closure problem: knowing S n (R) requires 
information about S n +\ etc. It is well known that arbi- 
trary closures of this hierarchy of equations are doomed, 
leading to predictions that are in contradiction with ex- 
periments. On the other hand no one succeeded to solve 
the hierarchy in its entirety. In contradistinction, the 
theory for the fully unfused F n does not suffer from this 
problem: it was shown that there exist homogenous equa- 
tions for F n in terms of F n , without any hierarchic con- 
nections to higher or lower order correlation functions. 
This fact allows one to proceed Jl^,|ll| to derive a variety 
of exact bridge relations between the scaling exponents 
of gradient fields and the scaling exponents of the struc- 
ture functions themselves, and to study the nature of the 
dissipative scales in turbulence, showing that in fact they 
are scaling functions with well defined scaling exponents. 

In considering the decorrelation times of many "fully 
unfused" space- time correlation functions we need to 
make a choice of which velocity field we take as our funda- 
mental field. The Eulerian velocity field won't do, simply 
because its decorrelation time is dominated by the sweep- 
ing of small scales by large scale flows. In || we showed 
that at least from the point of view of the perturbative 
theory one can get rid of the sweeping effect using the 
Belinicher-L'vov velocity fields whose decorrelation time 
is intrinsic to the scale of consideration. In terms of the 
Eulerian velocity ref. defined the field v(tq, to\r, t) as 

v(r ,t \r,t) = u[r + p h (r ,t \t),t] , (3) 
P L (n>,*o|*) = / u[ro + p L (r 0! to|r),r] . (4) 
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The observations of Belinicher and L'vov was that the 
variables v(r ,to\r,t) satisfy a Navier-Stokes-like equa- 
tion in the limit of incompressible fluid, and that their 
simultaneous correlators are identical to the simultane- 
ous correlators of u(r, t). 

Introduce now a difference of two (simultaneous) BL- 
velocities at points r and r' 

W(r , t \r, r', t) = v(r , t \r, t) - v(r , t \r', t) . (5) 

The equation of motion for W can be calculated starting 
from the Navier-Stokes equation for the Eulerian field, 

^+£(r, r', t)-i/(V? + V; 2 )] W(r , t \r, r', t) = 0. (6) 

We introduced an operator C(tq, to\r, r', t) 

C(r ,t \r,r',t) = PW(r ,to\r,r ,t) ■ V r 

+ P'W(r ,t \r',r ,t)-V' r , (7) 

where P is the usual transverse projection operator 
which is formally written as P = — V~ 2 V x Vx. The 
application of P to any given vector field a(r) is non 
local, and it has the form: 



[Pa{r)] a = J drP aP {i 



f)a (r). 



(8) 



The explicit form of the kernel can be found, for example, 
in ^ . In (0) P and P' are projection operators which act 
on fields that depend on the corresponding coordinates r 
and r' . The equation of motion (g|) form the basis of the 
following discussion of the time correlation functions. 

To simplify the appearance of the fully unfused, multi- 
time correlation function of BL-velocity differences we 
choose the economic notation Wj = W{ro, to\rjr'j, t): 

) = (Wi...W„) • (9) 

We begin the development with the simplest non- simul- 
taneous case in which there are two different times in (||) . 
Chose U — t + t for every i < p and U = t for every i > p. 
We will denote the correlation function with this choice 
of times as (r) , omitting for brevity the rest of the 

arguments. Compute the time derivative of ^F^\ with 
respect to r: 



p 

E 

3=1 



(Wi 



aw, 



dr g t 

3=1 

Using the equation of motion (|J) we find 



Wn> , (10) 



dr 



3=1 



(11) 

(12) 



with Cj = C(ro,to\rj,r'j,t). We remember that CjWj 
is a nonlocal object that is quadratic in BL-velocity dif- 
ferences, cf. Eq.(^). We reiterate that all the functions 
depend on 2n space coordinates that we do not display 
for notational economy. 

To understand the role of the various terms in Eq. jll| ) 
we will make use of the analysis of a similar equation 
that was presented in 1 1 1 1 . In that case we considered 
the simultaneous object F n = J-^ \{fS) and computed, as 
above, its i-time derivative. Obviously, in the stationary 
ensemble this derivative vanishes. Instead of \ (t) and 

Jn\( T ) w e got 2>2(t) and J ( ™[{t) which differ from 
the present ones only in that the summation go up to n 
instead of p. The crucial observations of |Tl[] are the fol- 
lowing ones: (i) ]im„_, JT„ \{t) = 0. For fully unfused si- 
multaneous correlation functions the viscous term in the 
balance equation disappears in the limit of vanishing vis- 
cosity, (ii) The integral in ™ \ (r) which originates from 
the projection operator converges in the infra-red and the 
ultraviolet regimes. This means that every term in the 
sum over j that contributes to 2? „ n \ (t) can be estimated 
as S n +i(R) j R when all the separation Rj = rj — r'^ are 
of the same order R. When we take the full sum up to n 
there exist internal cancellations between all these terms, 

(n) 

leading to the homogeneous equation D n [ = 0. 

In the present analysis the proof that the viscous term 
Sn\ ( T ) i s negligible when v — > is an immediate conse- 
quence of the previous result. Since we are taking only 
partial sums on j, the internal cancellation in2?£fj (t) dis- 
appears, and it has a finite limit when v — > 0. On the 
other hand JT„ i( T ) can only increase if we take r = 0. 

Thus again when v — > we can neglect «7^fi(T") in the 
fully unfused situation. Accordingly for v — > we have 



^ ) 1 (r)/dr + DW(r)=0 



(p), 



(14) 



The proof of convergence of the integrals in T>^ \ (t) fol- 
lows from the previous results, since the time correlation 
functions are bounded from above by the simultaneous 
ones. Accordingly, when all the coordinates are fully un- 
fused and the separations are all of the order of R, 



(p+i) , 



(15) 



Introduce now the typical decorrelation time t^\(R) 
that is associated with the one-time difference quantity 
3~n\( T ) when au the separations are of the order of R: 



(16) 



Remember that the simultaneous correlation functions 
of BL-velocity differences @ are identical || to the si- 
multaneous correlation functions of Eulerian velocity dif- 



J~n\( T ) = ^E( V J + V 1)< W i • • • W i • • • W "> » ( 13 ) ferences (H). Le - ^i P) i(°) = F n- Integrate (|14|) in the 



3=1 



interval (— oo,0), use the evaluation (Ua), and derive 
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RF r , 



JP+1) P 



(17) 



We see that from the point of view of scaling there is no 
p dependence in this equation: for different values of p 
only the coefficients can change. We thus estimate 



r„,i(i?) ~ RS n - 1 (R)/S n (R) oc R z 



(18) 



Here we introduced the dynamical scaling exponent z n _i 
that characterizes this time and found that 



ZnA 



1 + Cn-l — Cn ■ 



(19) 



Consider next the three-time quantity that is obtained 
from J- n by choosing U = t + n for i < p, ti = t + r 2 
for p < i < p + q, and £, = t for i > p + q. We denote 
this quantity as ( T i ' r 2)i omitting again the rest of 

the arguments. We define the decorrelation time rjf^ °f 
this quantity by 



J — oo 



^7(0,0). (20) 



One could think naively that the decorrelation time 
is of the same order as (|l8|). The calculation leads to a 
different result. To see this calculate the double deriva- 
tive of (ji.t-i) with respect to t\ and r 2 . This re- 
sults in a new balance equation. For the fully unfuscd 
situation, and in the limit c^Owe find 



&'^ ) (n,n)/dT 1 dv i + 2>^ 9 Vi,t 2 ) = , (21) 

where now 



v p+q 

< 2 9) (n,r 2 ) = ^ Yl 'Wi..X,W.....£,W/,, 



•Wn). 



In the RHS of @ we neglected two terms that vanish 

in the limit v — > 0. The expression for I?„ 2^ ( T ij r 2) con- 
tains two space-integrals that originate from the two pro- 
jection operators which are hidden in Cj and Using 
the same ideas when all the separations are of the order 
of R we can estimate with impunity 

< 2 9) (n,r 2 ) ~ J&Sf-ViTL,!*)/!? • (22) 

Integrating now Eq.(|2l]) over n and r 2 in the interval 
[— oo, 0] and remembering that 7"^^ (0j 0) = F n we find 



i?F n ~ -F 1 n+2[T n+2 2 



(p+l,9+l)i2 



(23) 



As before the scaling exponents are independent of p and 



q and we we can introduce a notation r„ 2 ~ r 



n,2 



Up 



to p, q-dependent coefficients 

[tuAR)? ~ R 2 S n - 2 (R)/S n (R) cx [R z ^] 2 . (24) 



We see that the naive expectation is not realized. The 
scaling exponent of the present time is different from (19): 



z n .2 = 1 + (Cn-2 - Cn)/2 



(25) 



The difference between the two scaling exponents z n ^\ — 
Zn,2 = Cn-l _ (Cn + Cn-2)/2. This difference is zero for 
linear scaling, meaning that in that case the naive ex- 
pectation that the time scales are identical is correct. 
On the other hand for the situation of multiscaling the 
Hoelder inequalities require the difference to be positive. 
Accordingly, for R L we have r„ 2 (i?) 3> T n ^(R). 

We can proceed with correlation functions that depend 
on m time differences. Omitting the upper indices which 
are irrelevant for the scaling exponents we denote the 
correlation function as T n ,m(j\ ■ ■ • T m ), and establish the 
exact scaling law for its decorrelation time. The defini- 
tion of the decorrelation time is 



(n . . . r m ) = [ (0...0). 



(26) 

Repeating the steps described above we find the dynami- 
cal scaling exponent that characterizes r„ !m when all the 
separations are of the order of R, r n rn oc R Zn < m : 



1 + (Cn- 



Cn)/m, , n — m < 2. 



(27) 



One can see, using the Hoelder inequalities, that z n ^ m 
is a nonincreasing function of m for fixed n, and in a 
multiscaling situation they are decreasing. The meaning 
is that the larger m is the longer is the decorrelation 
time of the corresponding m+l-time correlation function, 

Tn,p(R) > T n,q(R) for p < q . 

To gain further understanding of the properties of the 
time correlation functions, we consider higher order tem- 
poral moments of the two-time correlation functions: 



(28) 



The intuitive meaning of [r k )^\ is a k- order decorre- 
lation moment of J-^\(R, t), and note that its dimen- 
sion is (time) fe . The first order decorrelation moment is 

(v) 

the previously defined decorrelation time r„ j\ To find 
the scaling exponents of these quantities we start with 
Eq. (pT[) , multiply by r fc , and integrate over r in the in- 
terval (—oo,0). Using the evaluation (|l5| ) and assuming 
convergence of the integrals over r we derive 



t {p) k-i 

' n,l r 



-4 



f F^iydr, (29) 

J — oo 



where we have integrated by parts on the LHS. We stress 
that in deriving this equation we assert that k + 1 mo- 
ments exist; this is not known apriori. Using the defini- 
tion ( psl ) and for all the separations of the order of R we 
find the recurrence relation 



RS n (R)(r^l ~ S n+1 (R)(r^[ 



(p+i) 



(30) 
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The solution is 



i(p) 



(T n ,fe) 



R k S n -k(R) 

S n (R) 



(31) 



for fc < n — 2. The procedure does not yield information 
about higher k values. 

We learn from the analysis of the moments that there 

is no single typical time which characterizes the r depen- 

(%>) 

dence of There is no simple "dynamical scaling 

exponent" z that can be used to collapse the time de- 
pendence in the form \(r) ~ R^ n f(r/R z ). Even the 
two-time correlation function is not a scale invariant ob- 
ject. In this respect it is similar to the probability distri- 
bution function of the velocity differences across a scale 
R, for which the spectrum of ( n is a reflection of the lack 
of scale invariance. 

The main conclusion of this Letter is that all the dy- 
namical scaling exponents can be determined from the 
knowledge of the scaling exponents £„ of the standard 
structure functions S n (R)- All the scaling relations that 
were obtained above can be easily remembered using the 
following simple rule. 

To get the dynamical scaling exponent every integral 
over r in the definition of the decorrelation time ( f26[) , and 
every factor r in the definition of the moments (Eq) can 
be traded for a factor of R/W within the average of the 
correlation function involved. The dynamical exponent 
is determined by the resulting scaling exponents of the 
resulting simultaneous correlation function. 

The deep reason for this simple rule is the non- 
perturbative locality (convergence) of the integrals ap- 
pearing in £>„. Because of this locality one can estimate 
from the equation of motion (^|) 1/r ~ d/dr ~ W • V ~ 
W(R)/R. This means that we can use the substitutions 



t,J dr^ R/W(R) 



(32) 



as long as we use them within the average, and when all 
the separations are of the order of R. We propose to refer 
to this substitution rule as "weak dynamical similarity" , 
where "weak" stands for a reminder that the rules can 
be used only under the averaging procedure, and only 
for scaling purposes. The same property of locality of 
the interaction integrals was shown [|lO],[ll| to yield an- 
other set of bridge relations between scaling exponents 
of correlation functions of gradient fields and the scal- 
ing exponents £„. Those relations can be summarized by 
another substitution rule that we refer to as "weak dis- 
sipative similarity" ; It follows from equating the viscous 
and nonlinear terms in the equations of motion: 



W(R)/R. 



(33) 



Again "weak" refers to the reminder that we are only 
allowed to use these substitutions for scaling purposes 
within the average. Note that our weak dissipative sim- 
ilarity rule is weaker than the Kolmogorov refined simi- 
larity hypothesis which states the dynamical relationship 



(VW(i?)) 2 ~ W 3 /R. Both our rules are derived from 
first principles, while Kolmogorov's hypothesis is a guess. 

Finally, we should ask whether the results presented 
above are particular to the time correlation function of 
BL- velocity differences, or do they reflect intrinsic scaling 
properties that are shared by other dynamical presenta- 
tions like the standard Lagrangian velocity fields. The 
answer is that the results are general; all that we have 
used are the property of convergence of the interaction 
integrals, and the fact that the simultaneous correlation 
functions of the BL-fields are the same as those of the 
Eulerian velocities. These properties hold also for La- 
grangian velocities, and in fact for any sensible choice 
of velocity representation in which the sweeping effect is 
eliminated. Accordingly we state that the dynamical ex- 
ponent are invariant to the representation and in partic- 
ular will be the same for many-time correlation functions 
of Lagrangian velocity differences. 
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